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A cosmic axion, via the electromagnetic anomaly, induces an oscillating electric dipole for the 
electron of frequency m a and strength ^(few)xlO -32 e-cm, two orders of magnitude above the nu¬ 
cleon, and within a few orders of magnitude of the present standard model constant limit. We give a 
detailed study of this phenomenon via the interaction of the cosmic axion, through the electromag¬ 
netic anomaly, with particular emphasis on the decoupling limit of the axion, dta(t) oc m a —> 0. The 
analysis is subtle, and we find the general form of the action involves a local contact interaction and 
a nonlocal contribution, analogous to the “transverse current” in QED, that enforces the decoupling 
limit. We carefully derive the effective action in the Pauli-Schroedinger non-relativistic formalism, 
and in Georgi’s heavy quark formalism adapted to the “heavy electron” (m e >> m a ). We compute 
the electric dipole radiation emitted by free electrons, magnets and currents, immersed in the cosmic 
axion field, and discuss experimental configurations that may yield a detectable signal. 
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I. INTRODUCTION 

In the present paper we give a detailed study of the in¬ 
teraction of the cosmic axion through the electromagnetic 
anomaly, with the magnetic dipole field of an electron. 
For an oscillating cosmic axion field, we show that the 
electron acquires an effective oscillating electric dipole 
moment (OEDM). Our detailed analysis is subtle, and 
amplifies the result of our previous note [lj. The anal¬ 
ysis is perturbative, and is treated both classically and 
quantum mechanically. Not surprisingly, we find that the 
OEDM displays subtleties that are shared with the axial 
anomaly itself. 

We’ll assume that the axion fills the vacuum as a clas¬ 
sical coherent field, oscillating in a given frame with fre¬ 
quency m a , and may be associated with dark matter, as 
per the model in [2|. We will be interested in the ef¬ 
fect upon magnetic objects that are essentially at rest 
relative to the axion cosmic rest-frame. The axion field 
is designated as Oft) = a(t)/f a where a is the canoni¬ 
cally normalized axion field and f a the decay constant. 
The cosmic axion field is then 6(t) = 0ocos(ro a t) where 
m a ~ rrinf-nlf a is the axion mass. 

A simple hand-waving argument can be given for the 
existence of induced OEDM’s arising from magnetic mo¬ 
ments immersed in a cosmic axion field. Witten showed 
that a 0-angle in QED will cause magnetic monopoles to 
acquire electric charges proportional to 0, i.e., they be¬ 
come “dyons” a a. If we thus consider a pair of dyon 
and anti-dyon, separated in space by a very small dis¬ 
tance, we will have both a magnetic and an electric dipole 
moment where the electric dipole moment is equal to the 
magnetic moment times 0. At very large distances we 
cannot, for all practical purposes, discern the presence 
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or otherwise of the underlying magnetic monopoles, but 
the electric and magnetic dipole fields persist. A cosmic 
axion filling space is essentially an oscillating 0-angle, and 
we might expect by this argument, therefore, an OEDM 
proportional to Oft). 

An OEDM is indeed induced for the electron by the 
Feynman diagram of Fig.(l) where the solid dot vertex 
is the anomalous coupling of the axion to electromag¬ 
netic fields. Our result can be written as an effective 
interaction for the non-relativistic electron in the zero 
electron-recoil limit, with a time dependent electric field 
( e.g ., radiation or a cavity mode) as: 

/ d A x 2g a 0(t) pBohr X f yX(aO • E(x,t) (1) 

where g a is the axion-7-7 anomaly coefficient (defined 
below), pBohr the Bohr magneton, and ~^{x,t) is an ex¬ 
ternal radiative electric field. The induced effective os¬ 
cillating electric dipole moment is proportional to the 
magnetic moment, d e « 2g a OoCOs(m a t)piBohr ~ 3.2 x 
10 -32 (g a /10 -3 ) cos(m a t) e-cm. We use the term, “effec¬ 
tive,” because this arises at the level of a one-particle re¬ 
ducible Feynman diagram. Note that Eq.(l) is the limit 
of the action in the case of source-free, time dependent 
(radiation) fields, and one cannot naively take the limit 
Oft) —> (constant) without including the nonlocal terms 
in the full action, as in Eqs.(29), (30). 

Previously, an OEDM has only been considered to be 
specific to baryons. It arises, not by the electromagnetic 
anomaly, but rather directly via the QCD-induced ax¬ 
ion potential. The magnetic moment of the electron is 
much larger than that of the nucleon, and hence the 
axion-induced oscillating electric dipole moment is al¬ 
most three orders of magnitude larger than that of the 
nucleon. Since the current best limit upon any DC el¬ 
ementary particle EDM is that of the electron, of order 
d e < 8.7 x 10~ 29 e-cm, 1, the electron may provide a 
promising place to search for an oscillating EDM. 
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While the OEDM appears above to be proportional 
to 6it), there is, however, a catch: In the limit that 
dt9{t ) —> 0 (decoupling limit) the perturbative Feynman 
diagrams involving the anomaly must vanish. But how 
do we reconcile decoupling, i.e., derivative coupling of 
the axion, from a hard dependence upon 6(t)? 

The decoupling of the axion at zero mass is subtle. 
We will see presently that a nontrivial nonlocal term is 
generated by Fig.(l) that enforces the decoupling. This 
nonlocality is remniscent of the “transverse current” that 
arises in radiation gauge quantization of QED (see, e.g., 
section 6.3 of ref.Jg). The nonlocal term insures that 
the action, S{9, F ^ v ,...) can be brought to the form 
S'id^O, A v , ...)+(total divergence), where F 
This is a subtle property shared with the anomaly it¬ 
self whereby the manifestly gauge invariant form of S 
does not display the axion derivative coupling. However, 
upon integration by parts we can display the derivative 
coupling of the axion, as in S ', while relinquishing mani¬ 
fest gauge invariance. Since these actions are equivalent 
up to a total divergence, both the shift symmetry of the 
axion and the gauge invariance of QED remain valid in 
perturbation theory. Displaying the action as in S and 
taking the zero recoil limit of the electron, which is kine¬ 
matically valid in the m e » m a limit, we obtain eq. ©■ 
The price we pay for this symmetry is the nonlocality of 
the effective EDM action of the electron. The structure 
of the action is, however, determined completely by this 
symmetry, as we discuss in Section(III). 

Indeed, nonlocality arises even in the familiar case of 
a classical axion induced RF cavity mode. There, the in¬ 
duced electric field in the cavity, if ( t) satisfies a similar 
condition, if (i) = if(io) + c(0(t) — 9(to)). This happens 
simply because if(t) is governed by a first order inho¬ 
mogeneous differential Maxwell equation and requires a 
boundary condition. As we’ll see, the particular solu¬ 
tion of Maxwell’s equations for an induced electric field 
in a static background magnetic field is of fundamental 
importance in axion electrodynamics and drives most of 
the interesting phenomena. Modulo this subtlety, the ex¬ 
plicit calculation of the Feynman diagram as in Fig.(l) 
nontrivially confirms the argument based upon Witten’s 
dyons. The full calculations simultaneously provides con¬ 
sistency with decoupling via the nonlocal term. 

We begin by giving a detailed derivation of the ef¬ 
fective action of the electron OEDM in Section(II). We 
consider both the non-relativistic Pauli-Schroedinger for¬ 
malism for a resting electron, and also Georgi’s covari¬ 
ant heavy quark formalism for electron of 4-velocity v M 
[jlj. The latter formalism is adapted to the electron, 
which may be viewed as ultra-heavy in comparison to 
m a , and shows that the resulting interaction is of the 
form oc 9{t)i\> v o for t\) v = (1 + ^)^/2, with 
4-velocity v The results are consistent, and reveal the 
full effective action with the nonlocal term. 

In Section III.A we show that the structure of the ac¬ 


tion with the nonlocal term is completely determined 
by the axion decoupling, i.e., by the shift symmetry, 
a/f a —> a/f a + <j>, which is maintained in perturbation 
theory. While the physical effective action of the OEDM 
is consistent with the aj f a — > a/f a + 4> symmetry, we 
emphasize that there are no additional suppressions in¬ 
volving higher powers of m a , i.e., our OEDM physics is 
on par with the induced oscillating electric field in an RF 
cavity experiment. In Section III.(B,C), we observe how 
this nonlocality arises in well-known solutions to, e.g., 
the RF cavity experiments. 

To further probe this phenomenon, we show explicitly 
in Section IV that the classical Maxwell equations for 
a localized magnetic dipole, such as an electron in free 
space, leads to the emission of electric dipole radiation, 
i.e., the classical radiation field from a stationary elec¬ 
tron is that of a Hertzian electric dipole radiator. The 
classical calculation is compared to the quantum calcula¬ 
tion, and they are found to be consistent. Large magnetic 
fields imbedded in conductors likewise provide a source 
for such axion induced electric dipole radiation. 

In Section V we consider a possible experimental con¬ 
figuration for detection of this radiation based upon an 
array. This is a broadband simple radiator, and can be 
viewed as an array of high field magnets, or as a planar 
slab of conductor with a large magnetic field imbedded 
in the plane of the conductor. This can produce power 
output of upwards of order ~ 10~ 24 watts and appears 
to be detectable radiometrically. The main advantage 
over RF cavity experiments is that broad-band radiators 
do not require resonant tuning. We are encouraged 
by the simple estimates of the signal integration that we 
provide that this may lead to detactability, even in the 
challenging range 10 12 > f a > 10 10 GeV, or short axion 
wavelength. We note that there are several papers that 
touch on these and related ideas, e.g., i, 0 , and 0- 

Let us recall some basic concepts. The axion is a 
hypothetical, low-mass pseudo-Nambu-Goldstone boson 
(PNGB) that offers a solution to the strong CP prob¬ 
lem of the standard model, and simultaneously provides 
a compelling dark matter candidate. The expected mass 
scale of the axion is m a ~ m 2 // a where typical expected 
values of the decay constant f a range from ~ 10 10 GeV 
upwards 0 - 0 . 

The axion is expected to have an anomalous coupling 
to the electromagnetic field if- if, taking the form: 

f iu) F ^ F>iU = ~ 9a I (2) 

where F^ v = (1/2 )e l _ lvpa F p,J , and g a is the dimensionless 
anomaly coefficient. In various models we have 0 - 0 ]: 

g a « 8.3 x 1(T 4 DFSZ 

g a « —2.3 x 10" 3 KSVZ (3) 

In making quantitative estimates in Section V. we will 
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FIG. 1: Feynman diagram for axion induced electric dipole 
moment. Solid dot is the axion anomaly interaction, d~^ • ll 
The dashed line is the incoming axion field, 9. The solid line is 
the electron, where the electron-photon vertex is the magnetic 
dipole moment operator of the electron. 


use the KSVZ result (see Section V, in particular Table 
I, for a listing of our preferred parameters). 

Most strategies for detecting the cosmic axion ex¬ 
ploit the electromagnetic anomaly EMI together with 
the assumption of a coherent galactic dark-matter back¬ 
ground field Q, a/f a = 9(t) = 9 0 cos(m a t). In typi¬ 
cal RF cavity experiments such as ADMX, one applies 
a large external constant magnetic field to the cavity, 
if o (it suffices to apply this field only in the conducting 
walls of the cavity), and the anomalous coupling to 6(t) 
induces an oscillating electromagnetic response field, E r 
and ff r . The “cavity modes” can become excited, which 
can generate a resonant signal in the cavity. This offers 
the possibility of both detecting the existence of the ax¬ 
ion and simultaneously establishing that it is a significant 
component of dark-matter. We briefly discuss, for sake 
of comparison, the energetics of RF in Section III.C, and 
quantitatively in Section V. 

Recently several authors have considered alternative 
modes of axion detection [ 20 | - [ 22 | . in particular, the pos¬ 
sibility of observing an OEDM for the nucleons. In¬ 
deed, a small oscillating electric dipole moment for the 
nucleons is predicted with a frequency m a given by 
d/v ~ 10~ w 9(t) ~ 3.67 x 10 - 35 cos(m a f) e-cm [2l|. Thus 
far, this effect has only been considered to be specific to 
baryons. It arises, not by the electromagnetic anomaly 
oc g a , but rather directly via the QCD-induced axion po¬ 
tential. 

The axion induced magnetic moment of the electron 
is about two orders of magnitude larger than that of the 
nucleon. Since the current best limit upon any DC el¬ 
ementary particle EDM is that of the electron, of order 
d e < 8.7 x 10 ~ 29 e-cm, i , the electron may be a promis¬ 
ing place to search for an oscillating EDM and the axion 
itself. 


II. FEYNMAN DIAGRAM ANALYSIS OF 
INDUCED OEDM OF THE ELECTRON 

A. Nonrelativistic Pauli-Schroedinger Action 

We compute the axion induced OEDM of the electron. 
We go to the electron rest frame which we assume to also 
be the rest frame of the cosmic axion field, i.e., the frame 
in which the axion field oscillates in time with no spa¬ 
tial dependence, a/f = 9 0 cos(m a t). We use the Pauli- 
Schroedinger formalism, then follow with an analysis in 
Georgi’s heavy fermion formalism. 

The Pauli-Schroedinger Lagrangian is, of course, the 
first order Dirac action in an expansion in l/m e , and 
takes the form: 

■ (i$ - e~Pj ■ (i$ - e'Tj if (4) 

where if is a two-component spinor. Eq.(j4]) contains the 
term: 

~2 ' 0 "^) ' ( e "^) ^ 13 

(5) 

having used Oi<jj = Sij + ie^k&k- This is the standard 
magnetic dipole interaction, —^if'^if • B with g = 2 . 

Consider the time-ordered product of i x the magnetic 
dipole action with ix the axion anomaly action: 

~{i) 2 T J dfx-^if* eijkdiAjakif 

x| ^ dfy9^ Mp d a A fi d 1 A p (6) 

Assuming 9 has only time dependence, and integrating 
the anomaly by parts in time we have 

iT *4mJ dix ^ e ii^kifdiAjJ d 4 y(d 0 9)eim n Aid m A n 

(7) 

Note that this forces us into radiation gauge, as a 
Coulomb term would be ~ 9ei mn di(fd m A n and upon in¬ 
tegration by parts we would have ~ —di9eim n <fdmA n , 
which vanishes by the hypothesis that the axion field de¬ 
pends only upon time. 

Contracting the vector potential from the magnetic 
dipole interaction with either vector potential in the 
anomaly yields: 

J dAy J ^( x ) te b feCr fcV’(^) 

xdiG(x - y)d 0 9(y°)ejmndmA n (y) (8) 

where we’ve included a 2x combinatorial factor. Here we 
use the covariant Feynman propagator of the photon, 

i ~ G ( x ~y) ^ -ibijG{x - y)- 

d 4 k e ik d x -y) 


G{x - y) 


(27r) 4 k 2 + ie 


(9) 
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since only the spatial terms are relevant and the gauge 
dependent terms, oc A, are seen to cancel owing to anti¬ 
symmetries (see ref. 0 )- Our sign conventions and nor¬ 
malization are consistent with ref.[23j], where we have 
d 2 G{x — y) = 5 i {x — y). 

Using, Cijktjmn — $im$kn T dindkm■ we can write; 

= J d A y j d A x ip^c7 k ip{x) 

x diG{x - y)d 0 9{y°) ( diA k {y) - d k Ai{y )) 

( 10 ) 

Note that the result of eq.([8|) is manifestly proportional 
to d 0 9 and therefore the axion is derivatively coupled up 
a total divergence. However, when displayed in this form 
we do not see the manifest gauge invariance. This issue 
is addressed in detail in Section III. 

We assume that the electron is in a common rest frame 
with the axion. To a good approximation the electron 
operator is static, i.e., do{ip^a k ip) = 0 + 0(l/m). since 
the electron mass is enormous compared to the relevant 
momentum exchange, m e » m a ■ The electron can 
therefore absorb this tiny momentum without apprecia¬ 
bly changing its energy. 

We can therefore integrate the time derivative by parts, 
and drop d y oG{x — y) terms. We can also integrate the 
spatial derivative d Vk by parts, and note the transla¬ 
tional invariance of the propagator implies, d,,G{x—y) = 
-d x G{x - y). This yields: 

= tm j dAy J d ^ X ^ ak ^^^G( x ~ y)6 i.d 0 A k (y)) 

J d4y J d4x ^ a ki’( x )9kG{x - y)9 (didoA^y)) 

( 11 ) 

where V 2 = dkd k ■ Using the definition of the electric 
field in radiation gauge, Ei = —doAj, the static electron 
field permits us to integrate over a: 0 , f dx° V 2 G(x —y) = 
5 3 {x — y). We obtain the effective action (we have re¬ 
moved the i): 


is derivatively coupled. Our interaction behaves like the 
anomaly itself in its display of manifest gauge invariance 
vs manifest shift invariance.. 


B. Georgi’s Heavy Quark Effective Theory Applied 
to the Electron 

The previous analysis relied upon the non-recoil ap¬ 
proximation of the electron, which is justified since 
m e » m a . This limit implies that the 4-velocity of 
the electron is approximately conserved in its interac¬ 
tions with the comoving axion field. A more “covariant” 
representation can be derived, however, by using Georgi’s 
heavy quark effective field theory formalism [7j. We ob¬ 
tain the same basic structure as in the Pauli-Schroedinger 
case, but formulae involving the electric dipole moment 
will become the more familiar relativistic forms. 

We begin by defining heavy electron Dirac fields with 
fixed 4-velocity, v p : 

-a i 3 = i’v ip -a ip = 'Pv ( 13 ) 

The Dirac magnetic moment operator then takes the 
form: 

- ^IPv^^lpvF^ = (l/5£appy'y P v' y ) 1p v d a A 13 

(14) 

where the axion field is 9 = 0q cos{m a v p x^). 

We again compute the time ordered product of i x the 
magnetic dipole action with i x the axion anomaly action: 

-(*) 2 T J d A xip v {x) {-ise a p pl Y'v 1 )ipv{.x)d a A 0 {x) 
x J d 4 y9{y)e r ' K ^d v A K {y)d ul A p {y) 

= J d A y J d A xtp v {x) {^e a p pl Yv 1 )ipv{.x) 

xd a G(x - y ) {d v 9{y)^ p d u A p {y)) (15) 


ieg a 
2 TO 


d A x ip^a k ip 9(t)E k (x) 


+ J d 4 y J d 4 x d k ip^a k ip{x) G(x - y)9{y°)d i E i {y) 

( 12 ) 


Note again that propagator gauge terms do not con¬ 
tribute owing the the e-symbols. We now relabel, per¬ 
mute, and use the identities: 

e a / 3 7 P e Q?? “ J ' 5 = - (g V pQ%g p + (~l) p (permutations)) (16) 


The result we have obtained contains the local contact 
interaction, which is a conventional electric dipole form 
proportional to 9{t). It also contains a nonlocal compo¬ 
nent. If we choose to immerse the electron in a source 
free field, such as an RF cavity mode or light, then 
diEi(y) = 0 and only the contact term remains. 

We have written this expression in a manifestly gauge 
invariant form, and therefore we display the axion field 
without a derivative. Integration by parts, as in the 
general argument of Section III, shows that the axion 


and: 

- i'lpy^CTapIpv = 1p v (j 5 7aV/3 ~ ^'YyVa) 1p v (17) 

Note that we can integrate by parts and rearrange, to 
obtain: 

= ~2rn J d4y J ^ X ^vl5(y pl ipv [d 2 G{x - y){d p 9)A 7 

+d?d 1 G(x - y) {{d/ 3 9)A P - {d p 9)Ap)\ (18) 
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This form displays the manifest derivative coupling of the 
axion, but not manifest gauge invariance. 

Alternatively, we can manipulate eq. m, using 
d 2 G{x) = S 4 (x), and multiplying by —i to obtain the 
action: 

= ~Tm j 5<7p ^ ^ dFP1 

~ 7 ^ J ' d 4 y j d 4 x d p (ip v {l5<7fn)il>v{x)) 

xG(x-y)d 0 [6(y)F^{y)} (19) 

Here we have displayed manifest gauge invariance, while 
the derivative coupling of the axion is not manifest. We 
emphasize that ea. tfliil) and ea. fiTil) are equivalent up to 
surface terms. Note that in ea. flTTJl) we see the famil¬ 
iar covariant EDM operators in the heavy quark fields: 
1/A (75o> 7 )^0) ■ 

To compare the Georgi calculation result to the Pauli- 
Schroedinger case we go to rest frame 1 —> 0, and care¬ 
fully note sign conventions in ref. [23|: 

~~im J (OF a0 ) 

~^~2rri J d4x ^ a ^^ E ^ ( 20 ) 

and: 

J d * V J dAxdP ^v (75^)^- 

G{x — y)d~ t [6(y)F pl {y)\ 

-t J d 4 y J d 4 x di (i/Saiip(x)) 

G(x-y)-(0d j E j (y)) (21) 

where the latter spinors (without subscript, v ) are two- 
component spinors. These expressions are equivalent to 
eq. m- The result of eq. m is somewhat more general 
than ea. (fl2ll . the latter refering to the frame with electron 
and axion both having 4-velocity v p . 

III. ANOMALIES, OEDM’S, AND AXION 
DECOUPLING 

A. The General Form of the Action 

Presently we give a general argument as to why the 
nonlocal term occurs, and we will derive the general 
structure of the action by simply demanding the shift- 
symmetry of the axion. 

The effect we are discussing arises perturbatively, from 
the interaction of the axion anomaly with a magnetic 
moment and therefore it must share subtleties with the 
anomaly itself. The anomalous interaction of the axion 
with electromagnetic fields is written in eq.@: 

( 22 ) 


where 0{x) = a(x)/f a ■ 

We can perform a transformation that shifts the axion 
as 9 —> 9 + (j) where <j> is an arbitrary angle. If <f> is con¬ 
stant the theory is invariant under this shift, since the 
anomaly only causes the action of ea. (l 22 l) to shift by a 
total divergence, which has no effect in perturbation the¬ 
ory (a total divergence, or “surface term,” is oc the total 
incoming momentum less the total outgoing momentum 
of the Feynman diagram). Therefore, we conclude, a/ f a 
must be derivatively coupled in any perturbative process. 

We see this in eq. d^l) if we integrate by parts and 
write: 

- g a e^ P a j d 4 x d p ’6{t)A v F pa (23) 

where A M is the vector potential, and F^ = d^A^—d^Ap. 
Ea. (E51) displays the derivative coupling of the axion, 
however the manifest gauge invariance of electromag¬ 
netism is lost. Alternatively, ea. (l22ll maintains the gauge 
invariance of QED, but does not display the U(1)pq 
“shift” symmetry of the pNGB. Both symmetries are 
present in perturbation theory since ea. (l22l) and ea. (l23l) 
differ only by surface terms. We observed in Section II 
that this feature is shared by the Feynman diagram that 
generates a perturbatively induced OEDM. 

We generally define the covariant OEDM of the elec¬ 
tron (or similarly for any other object) in an arbitrary 
frame in which the cosmic axion field has four-velocity 
Ufj,, i.e., 6(x) = 6 q cos(u At a: A ‘), as: 

S' = g' J d 4 x e{x)S liV F pv (24) 

where is an antisymmetric odd parity dipole density 
(S°„ ~ ipo-fiu 7 5 , i/’; S pL , will be defined in terms of S 
momentarily). 

Since this result arises perturbatively from the 
anomaly, it must have the shift symmetry in common 
with the anomaly. In particular we must be able to move 
a derivative exclusively onto 9 via integration by parts: 

S' = - 2 g' J d 4 x d p 9{x)S^A v (25) 

We see that the gauge field, A p , now appears explicitly, 
exactly as happens in the case of the anomaly itself. 

However, in order for ea. (1^5l) to be valid, a constraint 
must be satisfied: 

d p S^ = 0 (26) 

A general solution to this constraint can be written as a 
nonlocal form: 

Sfiv = S°„(x) + J d 4 y d [p G(x - y)d x S° v]x (y) (27) 

(note antisymmetrization in y and v). Ea. (l26l) satisfies 
= 0 for any antisymmetric where the Green’s 
function satisfies: 

d 2 G(x — y) = 5 4 {x — y) 


9a 

2 e P vpi j 


d 4 x 0(x)F pv F pa 


(28) 
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and G(x — y) is as defined in eq.©. 

The action with eg. (1271) thus becomes; 


S' = g' J d 4 xd 4 yO(x)F pv (x) 

x S° tlv {x)8 i {x-y) + d [p G(x-y)d x S° ]x (y) (29) 


and S pl/ = ipa^'y 5 ^ is local. We can transpose the inte¬ 
grand, using d x G(x — y) = —d y G(x — y) and performing 
integrations by parts in y, to obtain: 


S' = g' J d 4 x 9(x)S° u {x)Ff MV 
+2 J d 4 x J d 4 y S° pu (x)G(x - y)d p d lp (9(y)F ll v] (y))( 30) 


Therefore, S' becomes a total divergence in the limit 
d t 9(t) —>■ 0. In particular, the first term is a total di¬ 
vergence in space for a spatially constant (9(f). 

In the case of a source-free electric field, V ■ E = 0, this 
becomes: 


g' J d 4 x 9(t)ft ■ ft (36) 

indistinguishable from a simple electric dipole moment 
interaction. Because ft = —dtA , we can integrate this 
expression by parts to write: 

g' f d 4 x d t 9(t)ft ■ ft (37) 


This is seen to be equivalent, e.g., to the Georgi formal¬ 
ism result of ea. (fl9)) after an integration by parts, where 
we have g' = eg a /4m and S 3 V = ip v 'Y 5 a fiu il) v . The role 
of the nonlocal term is therefore to maintain the shift 
symmetry of the axion. 

The nonlocal term reduces further in certian limits. In 
the limit of a constant 9(x) we have d p (9(y)F p u (y)) = 
Od^F^(y) = j v (y), where j v (y) is a source current for 
the electromagnetic field. Since d 2 G = S, the vector 
potential is given by: 


A-n{x) = J d 4 y G(x - y)j y (y) (31) 

Hence, for constant 9: 


S' 


9S^F^ - 


9S^F^ = 0 


(32) 


which is the simplest statement of decoupling. 

We can also consider a static limit for the electron in 
which: 


dy° G(x - y) = 


d 3 k e -^-^-v) 


J (2tt) 3 
1 1 
4-7T \x — y\ 


k 2 

= -i- (33) 

V 2 


where V 2 (l/i?) = — 47r<5 3 (i?) (e.g, see [gj] section 6.4). 
We find by explicit calculation of Fig.(l) that S pL , is the 
a local electric dipole density operator, oc a pv il)). 

For a purely spatially constant but time dependent ax¬ 
ion field, 9(t) and a nonrelativistic, static electric dipole 
moment, g'S^ = —g'S° 0 = | Pi(x ), where d t Pi(x) = 0 
this becomes: 

S' = g' J d 4 x 0(f) (ft ■ ft + ft ■ ft ft ■ ft^j (34) 
In an arbitrary gauge, ft = ft if — d t ft , we see that: 


S' = g' J d 4 x 6(t)ft ■ (ftp) 



and manifest gauge invariance is lost, just as in the case 
of the anomaly when the derivative is placed on the axion 
field. Eas. (l36l) and (l37l) differ only by surface terms that 
are irrelevant to perturbation theory. 

Note that if the electric field has a static source, such 
as an electric charge {e.g., an atomic nucleus), located 
at xo, ft -ft(y) = QS 3 {y — x) and, ignoring the explicit 
dt9 ~ m a terms in the nonlocal component of ea. dHOl) . 
then quasistatic action reduces to: 

g' J d 4 x (9(t) — 9{t — \x — xo\)ft(x) ■ ft(x) (38) 

Here we see that the nonlocal term subtracts a retarded 
axion field from the local value of the axion field. This 
difference vanishes as m a — 0. This is the way the de¬ 
coupling is generally maintained in solutions to Maxwell’s 
equations in classical configurations, such as RF cavities, 
as we see in the next few sections. This result, how¬ 
ever, implies that the electron OEDM may be difficult 
to probe in atomic experiments (such as the ACME ex¬ 
periment, 0) where \x — Xo\ ~ rsohr is the Bohr radius, 
since ( 9{t ) - 9{t —\x — f 0 |) ~ rn a r Bo hr « 1- 

We emphasize that a nonlocal operator structure in 
electrodynamics is not novel. It is encountered in the 
“transverse electromagnetic current” in QED, e.g., when 
we quantize in radiation gauge (see, eg, Q; section 6.3 
and eq.(6.28)). In that case the nonlocal term is essen¬ 
tial to maintain the causality of the theory in this gauge. 
The transverse current occurs when we have Coulombic 
sources and a nonzero, time dependent component of the 
vector potantial, A$. Aq has no time derivatives in the 
action and is therefore an instantaneously propagating 
field, and cannot represent a physical out-going on-shell 
photon. The equation of motion for Aq is V 2 Ao = —p(x), 
where p{x ) is a charge density. If we want to allow time 
dependent Aq, then V 2 c>oAo = —dop{x,t), but from cur¬ 
rent conservation we have dop = V • j where j is the 
3-current. Hence, we have doA 0 = —(1/V 2 )V ■ j. This 
means that if Aq is to be time dependent, then there 
must necessarily be a 3-current, hence a vector poten¬ 
tial, A. We impose the condition V • A = 0. A satis¬ 
fies (d 2 — V 2 )A — V(9qAo = j (the equation of motion 
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of Aq is unmodified by this). This is often written as 
(9 q — V 2 )A = jr- where jr is the “transverse current” 
0 which takes the form jr = j — V(1/V 2 )V • j. jx is 
identically conserved with the nonlocal term, is consis¬ 
tent with the radiation gauge condition V ■ A = 0, and 
thus gauge invariance is maintained. From this, Lorentz 
invariance is also maintained. 

Thus, introducing A 0 time dependence requires a non¬ 
local correction to the current to maintain a conserved 
current, hence gauge invariance. This is formally simi¬ 
lar to the structure seen in ea. (1271) which maintains the 
shift symmetry (analogue of gauge symmetry) of the ax- 
ion. The apparent nonlocality is arising because we are 
treating the “vacuum” as effectively containing a space- 
time dependence, through the background cosmic classi¬ 
cal 9{t) field. Physical amplitudes thus inherit a nonlocal 
dependence upon the history of the vacuum. 


B. Axion in an Infinite Volume Static Magnetic 
Field and Inherent Nonlocality 

The “integral form of axion decoupling,” as we have 
seen above arising from the nonlocal term, is a general 
feature of the solutions to the Maxwell equations in var¬ 
ious practical situations. Suppose we have an infinite 
universe which contains a uniform static magnetic field 
fio and a background oscillating classical axion field, 
a(t)/f = 9 (t) = 9q cos (m a t) (this can be considered as 
an infinite volume limit of an RF cavity experiment as we 
do below). We consider the sources for fi 0 to be far away 
from the region of interest and therefore fi xfi o = 0. 
The decoupling, m a —> 0 limit, becomes somewhat sub¬ 
tle, even in this case. We can analyze this classically. 

The axion anomaly will generate an electromagnetic 
field of the form: fi = fi r and fi = fi o + fi r where 
fi r and fi r are oscillating “response fields.” Maxwell’s 
Equations in these fields become: 

Maxwell (1): 


Maxwell (2) 


fi xfi r - d t fi r = ~g a fi 0 (8,9) (39) 


fi X fi r + d t fi r = 0 (40) 


and fi • fi r = fi ■ fi r = 0. 

The Maxwell equations are coupled first order inho¬ 
mogeneous differential equations. They are consistent 
with the decoupling of the axion as m a —> 0, since the 
source term is proportional to dt9. The vector potential 
in Coulomb gauge likewise satisfies: dffi r — fi 2 fi r = 
-g a % ( dt9) where fi r = —dtfi r and V • fi r = 0, a 
single second order inhomogeneous differential equation. 
For an infinite universe filled with the magnetic field fi o 
we have translational invariance in space. 


The Maxwell equations have a particular solution that 
is consistent with the symmetry of spatial translational 
invariance: 

fi r = g a ~fio [ d,Td T 9(r) = gafi 0 (9(t) - 9(t 0 )) 

Jo 

fir = 0 (41) 

The solutions necessarily require the specification of a 
single boundary condition at an inital time, which we 
take to be fi r (to) = 0 and fi r (to) = 0 . We em¬ 
phasize that this is a non-propagating solution (since 

fi r = 0) and represents a time dependent “dual rota¬ 
tion” O ifi 0 ^ fir 0- 

The dependence upon the initial condition introduces 
an apparent nonlocality, or “history”, into the observed 
fi r (t). Clearly fi x fi r = 0 hence B r = 0. Note that 
the vector potential is nonlocal, satisfying 

fi r =g a fi o [ dr' f d,Td T 0(T) (42) 

Jo Jo 

Our infinite universe with a static magnetic field and an 
oscillating axion field has acquired an oscillating electric 
field. This is an example of a general “theorem:” Any 
magnetic moment becomes an oscillating electric moment 
in the presence of the oscillating axion. The magnetic 
field need not be restricted to a constant all-space filling 
form as in this toy universe example; it can be, e.g., the 
local field surrounding a magnetic dipole moment, as we 
can see classically. The effect of the axion is to produce 
a time dependent electric dipole for the electron. 


C. Axion Induced Electric Field in an RF Cavity 

To detect a cosmic axion signal we can deploy a very 
large constant, externally applied magnetic field fig 
within a resonant cavity. The solutions to the Maxwell 
equations for the response fields will always involve the 
same particular solution we just encountered in the toy 
universe, but also now includes homogeneous solutions 
that are required to implement the boundary conditions 
of the cavity. 

Maxwell’s equations for the reponse fields are as in 
eas. (1391401) . We now have conducting boundary condi¬ 
tions at the cavity wall, r = R: fi(r = R) = 0 (and since 
the form of the S-field is parallel to the cavity axis, we 
have no constraint upon B(r = R)). 

With cylindrical coordinates, (p, <p,z), we find a cylin- 
drically symmetric solution: 

fi = (kJ 0 {pm a ) + g a ) B 0 9z 

fi = kJi(pm a )B 0 — (f (43) 

m a 

The electric field has the form of our free space solution 
proportional to g a, plus a homogeneous cavity mode so¬ 
lution proportional to k. 




We now apply conducting boundary conditions at the 
cavity wall, E(p = R) = 0. We thus determine k : 


t=- g JiPP- 


t=~ 


9a 


\ Jo (Rm , 

Ji(pm a ) 
Jo(R 


1 B 0 z9 



(44) 


Note again that the solution vanishes as m a —> 0 since we 
define 9(t) = f* dTd T 0(r), where to is an earlier time at 

which 111 {to) = 0. The solutions is therefore intrinsically 
nonlocaHn time. Henceforth we will omit the understood 
tilde on 9 , i.e., 9 —> 9(t) — 9(to). 

This is an idealized solution with a perfect resonant 
behavior, i.e., for a special cavity satisfying Jo(Rm a ) = 
0 the solution has an apparent infinite amplitude. Of 
course, in reality the amplitude is damped by dissipation. 
This results in a finite Q value, and modifies the solution 
to: 


t = 
t = 


( J 0 (pm a ) 

- g “ (—— - 1 


Bqz9 


dt9 

m a 


where 


F = \j (Jo(Rm a )) 2 + 1/Q 2 


(45) 


(46) 


Note that the total oscillating field energy in the cav¬ 
ity is oc Q 2 . The finite Q arises if we include a resistive 
damping term in the Maxwell equations. In an idealized 
perfect cavity the oscillating fields are the electric field 
90° out of phase from the magnetic field, and hence a 
vanishing time averaged Poynting vector. There is, how¬ 
ever, with finite Q there is an induced, small 0(1/Q ), 
temporal phase shift, that we have not written, such that 
and are not exactly 90° out of phase. This allows 
the Poynting vector at the walls of the cavity to average 
to a nonzero result. Hence power is extracable at a rate 
oc Q. If we attempt to extract power at faster rate then 
Q will decrease since the dominant power loss mechanism 
becomes the radiative extraction itself. We will consider 
some quantitative aspects of the RF cavity solution in 
Section V. 


IV. ELECTRIC DIPOLE RADIATION FROM A 
STATIONARY ELECTRON 

A. Classical Calculation 

First we consider the electric dipole radiation from a 
classical magnetic moment immersed in the axion field. 
This calculation has validity for intense classical mag¬ 
netic sources. It cannot be adapted to the case of an elec¬ 
tron which requires the quantum calculation. Nonethe¬ 
less, it is instructive to compute and compare it with the 


quantum case in Section IV.B. We will see that the radi¬ 
ation is generated by the physical OEDM. The magnetic 
dipole field surrounding the source, though it appears 
on the rhs of Maxwell’s equations as a source term in 
the axion background, does not itself radiate. Instead, 
this is associated with the non-radiating particular solu¬ 
tion encountered in the previous section, and it allows 
implementation of various boundary conditions at short 
distance for the electric field. 

The standard Maxwell equations with axion anomaly 
source in the presence of an arbitrary static, local mag¬ 
netic field Bo(r) are: 

Maxwell (1): 

^ x ~3 r (r,t) - dt^ r (r,t) = -go^o(r)dt9(t) (47) 


Maxwell (2) 


^ x ~3 r (r,t) + <9 t if r (r, f) = 0 


(48) 

and ^ = 0. 

Presently o(r) = ^ x ~^o(r) will be the static field 
of a classical solenoidal magnet, centered at the origin, 
f = 0, where can be written formally in terms of 
a magnetic dipole source rti 0 = rtiS 3 ^), (see ref.Q, 
Chapter 5.6). We have formally: 




hence: 


and: 


=- 


t 0 = 


^ x rfro 

(49) 

1 nt X 

47t r 3 

(50) 


1 

47T 


8ir 


rr$5 3 (r) 


— 


3 J? (J^ • r?t) 


(51) 

(see the magnetostatics discussion in Jackson, eqs.(5.55- 
5.64) fgj; Jackson states that this is a purely classical 
construction and cannot be unambiguously applied to a 
quantum mechanical electron). 

A subtlety arises in the present case with the particular 
solution encountered in Section III. Consider a “source¬ 
less dipole field,” i.e., one in which there is no rti5 3 (r) 
term: 


Bo = — 


1 

4-7T 


£)(*_*£■*)) (52) 

We consider Maxwell’s equations, replacing i^ 0 by B 0 on 
the rhs of ea. (H71) . We then have the particular solution 
to the vacuum Maxwell equations: 


and: 


t r (r,t) = - g a B 0 (r)9(t ) 


X(r,i) = 0 


(53) 


(54) 
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This is a solution since: 


0 = ^x(l|-3^<^) 


(55) 


(note that any potential severe surface term singularity 
arising here, e.g ., such as in i(rj /r 5 ), will be oc Sij or 
oc 7'i r ; -, but contracted with e^-fc from the cross-product, 
and hence zero). However, since lf r (r, t) = 0, this is a 
nonpropagating solution. It exists for any background 
static magnetic field satisfying x lf 0 = 0 . We will 
require incorporating this non-propagating solution into 
our full radiation solution momentarily. 

The existence of this solution has an important impli¬ 
cation: The sourceless dipole field surrounding the ori¬ 
gin, does not lead to radiation. The radiation comes only 
when we have a nonzero magnetic field if o which curls 
around a nonzero dipole source term. This is analogous 
to the infinite universe with a constant magnetic field ver¬ 
sus the RF cavity: it is the conducting wall of the cavity 
that enforces the boundary condition that produces the 
resonant magnetic and electric fields. 

We can now solve the Maxwell equations using re¬ 
tarded Green’s functions for a vector potential, A r (r,t), 
describing the radiative oscillating response fields in 
Coulomb gauge. The analysis mostly follows the text¬ 
book derivations as in ;(|, Chapter 9, but involves the 
non-propagating solution described above. In what fol¬ 
lows, we will pass to complex notation where 9(t ) = 
9q exp (im a t), and the physical response fields will be the 
real parts of the complex lf r and if r . The radiated elec¬ 
tromagnetic fields lf r (r, t) and lf r (r, t) are obtained as 
follows: 


if r {r,t ) = — ^-g a 0(t)exp(-im a \^\) 

4-7T 

(. . . /nt r (nt • "rh 

• Ml - exp(zm a | r |) + im a r\ ~ 3 -—- 

" rrl "r^(rrt • 


(56) 


and: 


~3r(r,t) = -^—g a d t 9(t) exp(—fm a |T*|) 
4-7T 


nl x ( — + 


im a l^ 


(57) 


One can readily verify that eqs. m E3 satisfy the 
Maxwell equations eqs. m USD with the source term 
—g a B 0 (r)(d t 9(t)). Notice that the second term in the 
brackets [...] in ea. dSSl) is the non-propagating solution of 
ea.(l53l) 

Taking the near-zone limit yields: 

lf r (f,i) —> 0 


which vanishes due to cancellation with the particular 
solution, and: 




1 

47T 


g a im a 9(t) 



(59) 


From the near-zone limit we can confirm the Maxwell 
equations and read-off the source structure: 


V x 


t r -d t t r = -±-g a d t 0(t)(^-^lti6 3 (A>) 


Til (T^-Th 

+ 73 - 3 ^5 


(60) 


Thus we see that the propagating radiation is due to 
the physical OEDM source, which induces the curling 
magnetic field, and is not due to the dipole magnetic 
field surrounding the source. 

In the far-zone we have: 

1W-,<) -3 

47r \ r r A r J 

(61) 


if r (r, t) 


gaml 

47r 


9(t) exp(—im Q |"7^|) 


Tit 


x 



(62) 


These are seen to be formally equivalent to the electric 
dipole radiation fields , where our magnetic moment m 
replaces the electric moment p in Jackson, @ eq.(9.18). 
Hence, the source of the radiation is the axion induced 
OEDM. 

From this we can compute the cycle averaged Poynting 
vector, if =< lf r x lf r >: 


t 


1 

32t r 2 


r, 2 ^n 4 fl 2 

9a m a9 0 




(7lt-l») 2 

Y 3 


(63) 


Using rrt = gBohrS , the angular differential emitted 
power, P, for our classical electron is therefore given by 
the classical dipole pattern (@, eq.(9.23)): 

dP 1 

h = 32 ^3>a« ohr S in2 0 (64) 

The total the emitted power is then: 

Plot = T^SaTO^oMBohr- ( 65 ) 


(58) 
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B. Quantum Calculation 

The quantum calculation is straightforward, but as one 
would expect is also somewhat tricky. We will summarize 
it presently. From the Pauli-Schroedinger result we have 
an effective action for radiation given by eq. m-- 

= ^m J (66) 

where we can neglect the nonlocal term since diEi = 0. 

The coherent axion field 9{t) = + e -trn a t^ 

has both incoming, e~ lmat and outgoing e +lmat com¬ 
ponents, and we drop the outgoing component for the 
process of conversion of an axion into a photon. Like¬ 
wise, we replace Ei(x) = doAi{x) by an outgoing pho¬ 
ton ikoeie +lkot ~ lk ' x with polarization e where e ■ k = 0. 
The electron is assumed to be localized in space and we 
write for initial and final two- 

component spinors Xi and Xf- 

Consider an emission amplitude of a photon in the xz 
plane (which corresponds to the polar azimuthal angle, 
<p = 0; we’ll integrate over <j> subsequently). The pho¬ 
ton 3-momentum is k = (sin(0), 0, cos(0)). The S 3 (x) 
distribution implies that the 3-momentum of the photon 
is arbitrary, constrained only by the on-shell condition 
and, finally, by energy conservation ko = \k\ = m a . The 
photon can, in principle, have two independent polar¬ 
izations, satisfying e t ■ k = 0, which we can take to be 
e*i = (— cos(0), 0, sin(0)) or €2 = (0,1,0). 

Consider first the case of the initial electron with spin 
up transitioning to a final electron, also with spin up, 

hence Xi = Xf = ^ 0 )' ^e 11 ! we see that: 

ei • Xf°Xi = sin(0) and e 2 • Xf°Xi = 0 (67) 

Therefore, only a photon of polarization e) is emitted 
in this case. So the amplitude for spin-up to spin-up is 
therefore: 

An = ^da^Bohrdoko sin(0) x 27t<5(to 0 - ko) (68) 

The corresponding transition rate is: 

t, !/ n ^ f k$sm 2 (d)d 3 k , ^ 

r tt = j{9ad-Bohrd o) J ( 27r )3 2 fc 0 ( 27r< H m ° ~ k 0)) 

1 f Tfl 

= j{ga9Bohrm a Oo) 2 J 9( 2 s ™?{0)d,0d<l> 

= X7r-(9a9Bohrm a Oo) 2 (69) 

1Z7T 

The emitted power is therefore: 

Ptf = m a T n = — ( g a 9Bohrm 2 a e 0 ) 2 (70) 

JLZ7T 

This is identical to the classical case of ea. (l65l) . 


However, there is, for a free electron, the possibility of 
a spin flip. Consider the case of the initial electron with 
spin up transitioning to a final electron, with spin down, 

hence Xi — q ^j Xf = ^ i ^ Now we see that 

ei • Xf&Xi = 0 and e 2 • Xf°Xi = * (71) 

Therefore, only a photon of polarization e 2 is now emit¬ 
ted. So the amplitude for spin-up to spin-down is there¬ 
fore: 


2 

An = x9a^BohrOoko2Tv5(ma - ko) (72) 

The corresponding rate is: 

1 f Tfl 

r tt = j(9a9Bohrm a 9o) 2 J 2^)2 sin {9)d9d(j) 

= x^-{ga9Bohrm a 9o) 2 (73) 

lD7T 

The emitted power is therefore: 

Pn = «taT n = T;( 9a9Bohrm 2 a 9of (74) 

Hence, in free space an electron will radiate with a total 
power given by 


7 

Ptotal = Pf^ Pfi = ( 9a9Bohr' n 'l a 9o ) (75) 

For an electron constrained to remain spin-up the power 
is that of ea. (17>5l) . 

In the following we will be interested in polarized elec¬ 
trons, such as electrons in ferromagnets. Polarized elec¬ 
trons are generally sitting in a polarizing H-field, and 
there is therefore an energy cost in flipping the spin. 
Hence, we are justified in dropping the Pjq rate in esti¬ 
mates of aggregated electrons in magnetic materials emit¬ 
ting dipole radiation. The quantum calculation, where 
we restrict the final state to spin-up for an initial spin- 
up state, gives identically the same power rate as the 
classical calculation. 

It is interesting that the axion field can also cause elec¬ 
trons to absorb photons. We will not further consider this 
“axion induced cooling” effect, but perhaps it, too, has 
some interesting consequences. 
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V. SOME QUANTITIVE ESTIMATES 

Presently we give some quantitative estimates for pos¬ 
sible detection schemes. We will briefly review the RF 
cavity, but then focus exclusively on the possible detec¬ 
tion of the radiation from coherent assemblages of mag¬ 
nets or large scale magnetic fields. We plan a more de¬ 
tailed treatment elsewhere l24j. 

We begin by defining a useful scaling parameter: rj = 
(/o/10 12 GeV). Axion parameters and conversions to nat¬ 
ural units are as follows: 


decay constant 

fa 

r, = (fa/ 10 12 GeV) 

axion mass 

m a 

6.02 x 10 “ 15 77- 1 GeV 

axion wavelength X a 

2irh 
m n c 

20.6 rj cm 

axion frequency v a 

c/X A 

1.46 x 10 9 Hz 

cosmic amplitude 

0 o 

3.68 x 10 - 19 p ff 

anomaly coefficient 

9a 

—2.26 x 10 -3 assumed 

magneton BBohr CGS 

eh 

2 m e c 

83.591 GeV -1 

1 statvolt/cm (cgs) 


6.92 x 10" 2 ° GeV 2 

1 watt 


4.09 x 10 " 15 GeV 2 

1 esla 


6.92 x 10 ~ 16 GeV 2 


Table I. Axion parameters and conversions . 1 2 


A. RF Cavity Energetics 

Let us estimate the signal power of a resonant RF 
cavity exp eriment. [This follows an estimate of Aaron 
Chou [25j pertaining to the ADMX experiment]. The 
AD MX microwave cavity is cylindrical and has a length 
of L ~ 10 2 cm, and radius of the cavity bore of R ~ 25 
cm, therefore a cross sectional area of ir (25) 2 = 1.96 x 10 3 
cm 2 . The volume of the cavity bore is V = 7 tR 2 L 


1 We assume the galactic halo energy density p g = 
PgalaxyK 0.3 GeV/cm 3 ) ; the cosmic axion field amplitude, 
a(t) jf a = do cos (m^t + <f>), where 0 is arbitrary. Note that 0q is 
independent of f a E3- 

2 One must be careful, as usual, in the definition of the Bohr 

magneton. In quantum electrodynamics we typically use the SI 
system of units, in which e^j/An = a = 1/137, and field energy 
density [Poynting vector] is (i? 2 + i^ 2 )/2 [7^ x 7^]. Note that 
we define the Bohr magneton [i Bohr Si — C5//i/2meC in the SI 
system (see eq[5j). If we choose CGS units where — ol the 

energy density [Poynting vector] is (7? 2 + il 2 )/87r [Ix'i /4tt] 
and the Bohr magneton is now HBohr CGS — €.cGS^/‘^ rn eC. 
The definitions in SI vs MKS differ by a familiar factor of 
VItt, i.e., PBohrSl = V^TTPBohr CGS- Put another way, our 
above SI calculation of radiated power, eq. ESP , yielded P = 
(!27T )- 1 {p. B ohr SI-) 2 = (127r)- 1 (v / 47rA‘Bohr CGS---) 2 - Had we 
used CGS field normalizations and computed the Poynting vec¬ 
tor directly we would have directly obtained the formula P = 
(3) 1 (p-Bohr CGS---) 2 - 


~ 1.96 x 10 5 cm 3 . We will assume Q ~ 10 5 , and an 
applied constant external magnetic field Bq = 7 Tesla. 

We see from eq. (1451) that the oscillating signal fields in 
an RF cavity are of order: 

l^rl - l^rl ~ 9aB 0 6 0 Q « 5.8 x 10~ 12 (cgs) 

(76) 

The total signal energy in the cavity is therefore: 

~ ^ (g a B 0 d 0 ) 2 xQ 2 VC ~ 1.85X10" 19 (Q/10 5 ) 2 ergs 

(77) 

where C is a form factor parameterizing the shapes of 
the cavity modes which we take to be C ~ 0.7. 

A damped driven simple harmonic oscillator of natural 
frequency uiq, driving force Fq exp (iuit ), satisfies, 

d 2 (j) + rd t (j> + u}q4> = Fo exp (iut) (78) 

hence: 

(j) = Re(cj> o exp (iuit))-, fo = ^ „ (79) 

(uJq — U} Z ) + lUJL 

On resonance this will have an energy stored in a cycle 
of = l^o^o | 2 / 2 , and energy lost in a cycle (which must 
be replaced in a steady state by the driving term), of 

Siost = (^r) x \\(f) 0 u>o\ 2 . We define: 


Q 


= 2 n 


£ 0 


Siost 


WO 

r 


(80) 


This holds for an axion RF cavity with ui 0 = m a . 

Much of what limits Q in an RF cavity is resistive lost, 
but we can define an accessible signal power output as: 


P 0 = £ 



m a S 0 



(81) 


Using m a = 27r x 1.46(/ Q /10 12 ) GHz = 9.1 x 10 9 Hz, we 
have a signal power Pq = 8.48 x 10 _ 22 (2e) watts, where 
our numerical coefficient corresponds to an efficiency of 
e = 1 / 2 . 3 

How challenging is the extraction of the signal from an 
RF cavity? We think this is challenging. Note that the 
signal power can be computed from the Poynting vector 
flux through an effective “aperature,” i.e., a hole in the 
cavity wall of area A out from which a tranversely polar¬ 
ized signal can be extracted, e.g., a “bung hole” in the 
barrel. From Maxwell’s equations with an an ohmic cur¬ 
rent in the cavity wall, one finds electromagnetic fields 
that attentuate over a skin depth <5 ~ 1/ ^/m a a/2, and 


3 A. Chou [25| uses m a =27r x 1.0 x 10® = 6.28 x 10® GHz, 

Q = (1.3) x 10 s , and a magnetic field derived from the total 
magnetic field energy, Ub = 4 MJ, or Bo = 7.16 Tesla, yielding 
8.4 X 10 —22 watts, consistent with rescaling our above result. 
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which have a small non-zero Poynting vector that is re¬ 
duced by a factor of Q relative to the cavity energy den¬ 
sity. The radiated power out of the cavity is then: 


P' 

M) 


1 &lA ou t 

Q V 


(82) 


Given these two routes to computing the power, Bo an d 
Pq , we can compute the ratio of the aperature area, A out 
to the surface area, Aq, of the cavity in terms of e. We 
readily find: 


A 0 ut 

Aq 


e 


(83) 


If we choose e = 1/2, then A out /Ao ~ 0.5 which is 
large, and implies that extraction of energy through a 
physical aperature is limited, since most of the surface 
area of the cavity is in the aperature hole! We can 
extract signal through a non-perturbative, small aper¬ 
ature area, of order A out ~ 7t(A q /4) 2 , i.e., a “quarter- 
wave aperature.” For the Poynting vector magnitude 
eci. (15i?l) apropos ADMX, we then have Pq = 7T-R(A a /4) 2 = 
2.36 x 10~ 24 77 2 (Q/10 6 ) watts, and e = 4.7 x 10~ 3 . We 
would thus take a significant hit in the output power with 
a less perturbative aperature size. 

Hence, if appropriate impedance matching of the cav¬ 
ity to an output receiver can be acheived, with e ~ 1/2, 
one would be able to attain power output of the order of 
~ 10 -21 from an RF cavity. The remaining bottleneck is 
then probing for a signal, which requires integrating and 
analyzing the output and searching for a signal/noise ex¬ 
cess for a given physical cavity tuning. This can be done 
on the order of minutes, but then the cavity must be 
retuned to a different resonant frequency and the signal 
integration process repeated. 


The signal fields in an RF cavity are driven by the par¬ 
ticular solution to Maxwell’s equations with the axion, 
which induces an oscillating electric field, ~^ r = g a 0 o^o, 
and a vanishing magnetic field ~$ r = 0 throughout space, 
including within the conducting cavity walls. The partic¬ 
ular solution has vanishing Poynting vector and cannot 
propagate, since ~^ r = 0. The resonance is determined 
by the presence of the conducting wall of the cavity. In 
this wall the induced electric field g a 0 qI^q generates a 
physical current. This physical current, in turn, stimu¬ 
lates emission of radiation into the cavity, with the effect 
that oscillating electric and magnetic modes, l^' r and if/, 
are generated in the cavity. The net electric field at the 
wall vanishes, lf r (R) + E' r {R) = 0. At resonance the 
induced fields are oc Q and the signal energy density is 
oc Q 2 . Finite conductance in the wall causes a slight 
phase shift of if/ and if/ away from 90° by an amount 
oc l/Q, leading to a the non-zero time-averaged Poynting 
vector of order Q. This is what we depend upon for a 
detectable signal. 


Note that the bulk magnetic field, Bo, within the inter¬ 
nal volume of the cavity, generates the particular solution 
there, but this is playing no role in the physics — only 
the B o field in the walls of the cavity is relevant! It would 
be sufficient to have the magnetic field localized only on 
the wall and not filling the internal volume. Moreover, 
this means that if we can apply a large magnetic field in 
any conducting material (or dielectric) it will then radi¬ 
ate observable power into the vacuum. This is the basis 
of an array radiator we now discuss below. 


B. Electric Dipole Radiation from Magnets 

We now turn to the radiation produced by localized 
magnetic fields, i.e., induced oscillating dipole moments, 
in the cosmic axion field. Using the formula of ea. flBSl) . 
(where we include the factor of \/Att into pBohr, or, use 
the CGS formula, P = ( g a 9o PBoh r m 2 ) 2 /3 directly) we 
can estimate the power emitted by a single electron in 
free space, immersed in the cosmic axion field: 

P e ~ 5.192 x 1(T 82 (?7)- 4 watts. (84) 

This follows the dipole distribution of ea. (l64l) . 

While this is infinitesimal power, it will add coherently 
for an assemblage of many polarized electrons within the 
near-zone of the radiation field, e.g., in an approximately 
a quarter wavelength volume, ~ (A a /4) , which we will 
define to be a “unit cell.” The unit cell can be viewed as 
a Dirac 5-function source, and we can construct an array 
of such sources. 

For a mole of polarized electrons comprising a unit 
cell, within the near-zone, the power becomes enhanced 
to: P Fe = N 2 vogadro P e ~ 1.88 x 10 - 34 (t 7 )- 4 watts. More 
generally we can assemble many moles of polarized elec¬ 
trons within the full quarter-wavelength volume to maxi¬ 
mize output power. The quarter wavelength volume, de¬ 
fined by the axion wavelength, is Vo ~ (5.216 g cm) 3 = 
1.420 x 10 2 77 3 cm 3 . For iron (Fe), of atomic mass A = 56 
gm/mole and density p ~ 7 gm Fe/cm 3 , we can therefore 
assemble pVo/A = 17.74ry 3 moles within a quarter wave¬ 
length volume. Assuming one polarized electron per iron 
atom, the power emitted by such a magnet is then of or¬ 
der (pVo/A) 2 P = 5.93 x 10 — 32 ? 7 2 watts. Note the scaling 
behavior in g is now controlled by overall the power per 
electron, g~^, times the quarter wave volume squared, 
(g 3 ) 2 for a net scaling oc g 2 . 

Note that the product of the external dipole mag¬ 
netic field with approximate volume it occupies, -BoVo, 
can be used in place of Npsohr i n computing the 
radiation fields. The Poynting vector is then K = 
(127t) _1 (g a dom 2 a BoVo) 2 , valid up to a Vo of order a quar¬ 
ter wavelength. The results using this formula are com¬ 
parable to the coherent multi-electron calculation. 

Neodymium-ferrite magnets can produce fields up to ~ 
1.4 Tesla. A unit cell composed of this material (using 1.4 
Tesla xVq in place of fiBohrNAvogadr oN m in the estimate 
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for Fe) yields a power of about Pn<i ~ 1.89 x 10 -30 77 2 
watts. Clearly, there is a significant advantage to larger 
field strengths. 

Bear in mind that arrays of unit cells spaced external 
to the near-zone produce a Poynting vector flux that is 
enhanced by N 2 , but this is focused within a solid angle 
of ~ 1/AT, i.e., the net power output is not coherently 
enhanced as N 2 for larger arrays, but only as N. We 
can, nonetheless, use an array of N unit cells to enhance 
the overall power output of the array oc N and focus the 
energy into a small solid angle to a receiver. 

We can also consider a “unit slab cell.” This is plate of 
conductor lying in the xy plane with a strong magnetic 
field that is imbedded and polarized in the plane of the 
conductor, e.g., in the y direction. Such a plate will radi¬ 
ate coherently in the directions normal to the xy surface, 
in the z direction, and the radiation will be electrically 
polarized in the y direction. 

In Appendix A we derive the solution for the ra¬ 
diation field from maxwell’s equations. The radiation 
from the xy slab into the z zenith direction is I^ r = 
—gBodo cos (m a t — kz)y and ~§ r = gBodo cos (m a t — kz)x, 
and the time averaged Poynting vector is given by: 

P 0 = -L(gB 0 e 0 ) 2 Az (85) 

87T 

Here A is the surface area of the slab and a factor 
of 1/2 arises from time averaging the Poynting vector. 
Given that the thickness of the conducting surface ex¬ 
ceeds the skin depth, the upward going radiation matches 
the gBodo particular solution in the skin of the conduc¬ 
tor, which sources the emitted radiation. Note that here 
we have assumed a negligible B 0 above the plane of the 
radiator. If Bq is constant in z above the plane of the 
radiator, then the radiated power will acquire an inter¬ 
ference term and ea. (|85l) will receive an overall factor of 
(1 — cos (kz)). This oscillatory behavior in space could 
prove useful in diagnosing a signal. 

We will therefore define a “unit slab cell” as a 1 mx 1 m 
conducting plate ( e.g., copper) with a 1 Telsa field in the 
plane of conductor. While this may be an unrealistic 
construct experimentally, we can use the result to scale 
to other parameters. We thus obtain for the unit slab 
cell, P s iab = 8.27 x 10“ 29 (Bo/1 T) 2 ( Area/1 m 2 ) watts. 

Essentially what we are doing presently is “opening 
up” the RF cavity, and creating a simple radiator. We 
will not be storing energy in a cavity, hence Q = 1 for 
us. Instead, by constructing an array of slab cells we 
would be exposing the largest possible conducting surface 
area and coherently focusing the axion induced emitted 
radiation toward a receiver antenna, which lies at some 
point in the z direction. 

Such a plate may not be easy to construct, but we 
can scale in area and field strength to match engineering 
constraints. This could be constructed e.g., by placing 
race-track magnetic windings wrapped around the plate 
in the x direction; alternatively, the plate could be seg¬ 
mented into x direction strips spaced by the solenoids 


aligned in the y direction. Segmentation of this config¬ 
uration is tolerable provided the average large distance 
(compared to an axion wavelength) is the regular square. 
In summary: 


Single electron 
Fe Unit Cell 
1.4T Nd Unit Cell 
lm 2 IT slab cell 


P e ~ 5.192 x lO -82 ^) -4 watts. 
Ppe ~ 1.88 x 10 — 34 (? 7 ) —4 watts 
PncL ~ 1.89 x 10 -3 Q t 7 2 watts 
Psiab = 8.27 x 1(T 29 watts 

( 86 ) 


For magnets, we expect to encounter complications 
due to the conductivity of the material, and the effects 
of eddy currents that are generated there. Most high 
field magnets are good conductors at GHz frequencies, 
and one might expect that the electric field will be nul¬ 
lified by the response in the material and hence in the 
near-zone of the radiation field. However, if one inspects 
eq. (1561) one sees that the electric field is vanishing in the 
near-zone, due to a cancellation of the electric dipole ra¬ 
diation with the particular solution from the oscillating 
axion field throughout space (analogous to the cancella¬ 
tion in the wall of the RF cavity). The radiation field 
in the near-zone is therefore predominantly due to the 
curling magnetic component of eq. (EH EU). Therefore, 
the emitted power results for aggregate magnets of or¬ 
der a quarter wavelength in size is a collective effect and 
localized eddy currents may not be an appreciable effect. 

One might think that we could simply replace the 
cylindrical magnets with small solenoids. Here the prob¬ 
lem arises that the solenoidal H-field itself can only gen¬ 
erate the non-propagating particular solution in the vac¬ 
uum, and we require the interaction of this field in matter 
to produce the radiation (this is akin to the RF cavity 
case). The solenoidal magnet is clad with the conduct¬ 
ing wire windings. The magnetic field can only penetrate 
significantly into this material within the interior of the 
solenoid, to generate a small cavity radiation (far below 
resonance) within the magnet. The external field is weak 
at the conducting material surface and does not tend 
to lie within the skin depth of the windings there, and 
we do not expect significant radiation to be generated 
by a solenoid itself (this is equivalent to the fact that a 
sealed RF cavity tends not to radiate externally if the 
return flux at the external conducting surface is small). 
Nonetheless, we could exploit strong field solenoids by 
arraying them in the xy plane with a conducting mate¬ 
rial strategically spanning the inter-magnet space (e.g., 
“fins”) in the xy plane. This will produce a coherent ra¬ 
diation signal in the z direction, which we discuss below. 
The number of unit-cell magnets in the array N then 
becomes equivalent to the RF cavity Q. 

Given an xy array of magnets with a polarization in 
the x direction, we can compute the discrete sum to ob¬ 
tain the array radiation fields and Poynting vector. Al¬ 
ternatively, we can take a continuous limit and the array 
becomes equivalent to a slab of conductor with an imbed¬ 
ded magnetic field in the x direction. We can imagine 
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using race-track solenoids wound around the slab, with 
exposed segments of conductor as the radiation surfaces. 
Other configurations can be investigated. In the follow¬ 
ing we estimate the signal and integration time necessary 
to detect in this kind of general xy array configuration. 


C. Axion Signal in a 2-D Axion —>• Photon Antenna 

We presently consider a schematic experimental con¬ 
figuration. This is a variation on a scheme proposed by 
flOl ]. The purpose of this is to provide an example of 
how one might attempt to exploit axion induced electric 
dipole radiation for detection in a broadband antenna. 
This is a rough “initial pass,” and we will refine it else¬ 
where [24| . 

Consider a smooth surface (a floor) which we define to 
be the xy plane. In this plane we assume we have placed 
a number N of the unit cells defined in the preceding 
section. The power emitted into the z direction is then 
given by: 

Pslab = TT - (gaBo^o ) 2 Atotal (87) 

57r 

where A tota i = NA is the total exposed conducting sur¬ 
face (and, as discussed above, if B 0 is constant above 
the plane of the radiator the result becomes P s iab = 
-^(g a B o 0 o ) 2 Atotaiil - cos (kz)) ). 

Details of geometric focusing to a receiver antenna will 
be developed elsewhere. It is not hard to see that the 
beam of radiation from the xy array will be coherently 
focused into a small solid angle and can be collected by 
a parabolic antenna looking down on the array (A vari¬ 
ation on this 0 would be to arrange the slab cells on 
the surface of a mathematical sphere of radius R and col¬ 
lect the signal at the focal point of the sphere; we think it 
may be advantageous to use an independent parabolic an¬ 
tenna and a planar array). The issues of antenna design 
and optimization are beyond the scope of our present dis¬ 
cussion. Let us crudely estimate presently the expected 
signal sensitivity and baseline requirements. 

Quantitatively, for an extremely optimistic “bench¬ 
mark,” we’ll assume Bq = 10 Tesla and that we have 
configured an array of 10 x 10 slab cells into an effective 
10m x 10m = 10 6 cm 2 conducting surface. We therefore 
we have an emitted power of P ar ray ~ 8.27 x 10 - 25 watts 
that is independent of r). This is, as we’ll see momen¬ 
tarily, more than adequate for detection in a cryogenic 
environment, and we will illustrate the scaling with ar¬ 
ray area, temperature, and magnetic field to back down 
to a more minimal experimental scale. 

The signal received in our antenna will be electronically 
filtered into a “pass-band” that we take to be of order 
~ 0.5m a to ~ 1.5m a for a hypothetical m a = 1.46 x 
10 9 77 -1 Hz. The pass-band is subdivided into frequency 
bins A/. These bins are given by the natural linewidth of 
the axion itself, defined by the so-called “axion Q a . v The 
axion Q a ~ 10 6 arises as a line broadening due to motion 


within the galactic dark matter halo, which generates 
fluctuations in the axion velocity, (3 ~ 10~ 3 . Hence, our 
frequency bins are of order A/ ~ m a x /3 2 ~ 10 3 Hz, and 
we can select ~ 10 6 such bins within our pass-band. The 
pass-band signal is recorded by the radio receiver over a 
long integration time and Fast Fourier Transformed. The 
experiment can be repeated for other pass-bands (other 
hypothetical m a ’s). 

The thermal noise power in a bin is given by the (one¬ 
dimensional) Bose-Einstein distribution for thermal pho¬ 
tons incident on the the receiver. The noise power in 
a single frequency bin is P no ise = 2TA/. We will as¬ 
sume cooling of the antenna (and other noise sources) 
to T = 0.1 °K corresponding to P no ise = 4.02 x 10~ 21 
watts. Hence our signal/noise ratio is Parray/Pnoise = 
2.06 x 10 -4 . This is small, but the number of noise pho¬ 
tons random walks in the number of “samples,” tAf, 
hence a signal can be observed with sufficient integration 
time to reduce the fluctuations in the noise, “A N noise ” = 
(Pnoise/'m-a)\/ tAf, below the number of signal events, 
A signal = (P 0 /m a )tAf. 

This is summarized by the Dicke Radiometer Formula 
for t: 


t = 



( 88 ) 


If we plug in the results for the minimal array we obtain: 


t = 9 hours 


(89) 


We give a tabulation of integration times for various ar¬ 
rays, each assumed to be 10m x 10m. 


Array 

Bo 

power 

time 

ND-Fe 

1.4 Tesla 

7.6 x 10 -27 watts 

3.2 years 

Solenoids 

3 Tesla 

6.8 x 10 -26 watts 

56 days 

10m x 10m 

10 Tesla 

8.27 x 10~ 25 watts 

9 hours 


TABLE II. Signal integration times for planar arrays of 
conductors, 10m x 10?n. 


For example, we can use fixed field ND magnets, with 
magnets spaced a half-wavelengths to form a fixed field 
array. We also consider 3 T solenoids, forming a layer 
under a conducting plane, where the return flux lies in the 
plane (various other configurations can be considered). 
Note the scaling in each case is: 


t oc 


(BoV /10m x 10ro\ 2 / T 
V"B ) ^ A ) \0dAA 


2 


(90) 


For example, a 20m x 20 ?tt Nd-Fe magnet array we 
would have an integration time reduced by l/(4) 2 , or 
2.4 months. Such an array would require ~ 4 x 10 4 mag¬ 
nets of order ~ 5 cm in scale length spaced by ~ 10 cm. 
Note that, despite the choice of quarter wavelength unit 
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cells magnets and spacing of order A Q /2 for some 77 ~ 1, 
the signal in the z direction is fairly broad-band and in¬ 
dependent of 77 ~ 770 - In Appendix A we describe briefly 
how the discrete array power approaches the slab power 
in the continuum limit. 

Note that if the xy array is a slightly spherical surface, 
then dipole radiation will be focused to the focal point, lo¬ 
cated at the mathematical center of the sphere. The cells 
would be emitting dipole radiation to the focal point in a 
common phase (see [lfl])- The use of a parabolic reflector 
over a flab xy planar array has the equivalent common 
phase relationship from the a flat array at the receiver 
focal point. Nonetheless, one might exploit the combi¬ 
nation of spherical array and semi-parabolic receiver an¬ 
tenna. 

With an observed signal in hand, one would be able 
to exploit interference effects such as the oc (1 — cos (kz)) 
modulation above the array in aconstant Bq field. It 
may also be possible to exploit polarization of the radia¬ 
tion and multiple receivers to reduce signal to noise and 
cryogenic requirements. The signal will be electrically 
polarized with the oreintation of the magnets or current, 
while the noise is unpolarized. 

This is a simple radiating system that does not de¬ 
pend upon a resonance condition and does not require a 
physical retuning of the array, and is fairly broadband. 
A more detailed discusion of a conceptual array-based 
axion search experiment will be given elsewhere 1241 . 


VI. CONCLUSION 


Our principal result is that the perturbative interac¬ 
tion with the cosmic axion, though the anomaly, causes 
an electron to acquire an oscillating electric dipole mo¬ 
ment. The result is general: any static magnetic source 
field in the presence of the cosmic axion will generate 
an oscillating electric field, hence any magnetic moment 
becomes an effective electric moment 4 . 

As an explicit example, we demonstrate that classi¬ 
cally that a stationary magnetic dipole field radiates as 
an oscillating electric dipole field; the near-zone struc¬ 
ture is that of a time dependent “Hertzian” electric dipole 
source of frequency to q . Likewise, the electron will have 
an OEDM that is proportional to its spin, and will ex¬ 
perience electric dipole interactions with applied electric 
fields. In addition we have the appearance of the nonlocal 
term, the analogue of the transverse current in radiation 
gauge. This nonlocal effect is present in any axion elec¬ 
tromagnetic interaction, such as the case of fields induced 
by the axion and applied magnetic fields in RF cavities 


4 As noted in ref.flj, the converse is not true, 2 .e., a static electric 
field does not become an oscillating magnetic field, a consequence 
of the fact that the cosmic axion field effectively breaks Lorentz 
invariance. 


(see Section V.). It does not affect the oscillating elec¬ 
tric dipole interaction of the electron with a source free 
electric field, i.e.as in radiation or cavity modes. 

We have obtained an induced oscillating electric dipole 
moment for the electron, proportional to the magnetic 
moment, 2g a 9o cos(m a t)fJ.Bohr- The result is quantita¬ 
tively « 3.2 x 10 -32 x 10 - 32 (< 7 a / 10 -3 ) cos(m a t) e-cm. 
The result is two orders of magnitude greater than the 
typical result expected for the nucleon, dw ~ 3.67 x 
10 -35 cos (m a t) e-cm [ 2 l|, about four orders of magnitude 
from the limit on the constant electric dipole moment of 
the electron, d e < 8.7 x 10 ~ 29 e-cm, [||. 

The result is a general low energy theorem and applies 
to any static magnetic system. Axion electromagnetic 
anomaly effects are essentially local oscillating dual ro¬ 
tations that lead to potentially observable signals. The 
axion anomaly perturbs the system by locally producing 
a physical, infinitesimal, time dependent dual rotation of 
<51^ = g a 77 9o(t)B, effectively rotating a magnetic mo¬ 
ment to an oscillating electric dipole moment. The dual¬ 
ity of axion electrodynamics also implies the absence of 
induced oscillating magnetic moments from static electric 
dipoles. 

The requirement of an explicit appearance of dtO(t) in 
physical quantities is a source of confusion to some peo¬ 
ple. This issue is not exclusive to OEDM’s, but also arises 
at the classical level in well-known solutions to axion- 
Maxwell’s equations, such as in conventional RF cavities. 
As shown in Section V., in an RF cavity with an applied 
classical constant background field Bq , an oscillating elec¬ 
tric field develops that has the form E oc 9(t)Bo. One 
does not see the explicit dt9(t) in the physical E, so one 
might wonder how it turns off in the m a 0 limit? The 
Maxwell equations for E and B are of the linear, first 
order, inhomogeneous form, and must therefore have a 
boundary condition. If one is careful, therefore, one finds 
that the electric field is proportional to: 


6{t) = dr 9 t 0(t) = 6{t) — 9(to) 


(91) 


'to 


Hence, 6(t) 0 when dtd(t ) —> 0, owing to the boundary 

condition. So, when one measures E(t) in an RF cavity 
one is only measuring its value relative to an earlier value 
E(to). In the present case of OEDM’s, however, this 
nonlocality is more subtle, but in the low energy zero- 
electron recoil limit of interest it reduces back to oc 9(t). 

The existence of such phenomena may imply a number 
of potentially sensitive venues. The collective magneti¬ 
zation of any subtance, e.g., polarized ferromagnets such 
as iron or strong rare earth magnets such as neodymium- 
ferrites, will acquire induced oscillating electric dipole 
moments that may provide accessible signals. 

We have presented estimates for some xy array config¬ 
urations that probe the coherent radiation generated by 
an assemblage of a large number of polarized electrons 
or electruic currents. Certainly other geometries can 
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be considered, including ID and even 3 D “crystal” ar¬ 
rays. These configurations can be built out of “unit cells” 
which are quarter (axion) wavelength volumes composed 
of magnetic materials. We find that reasonable labora¬ 
tory scale configurations, particularly the 2-D large array, 
may provide detectable coherent power outputs in excess 
of ~ 10~ 25 watts. These arrays have signal/noise integra¬ 
tion times that shrink inversely with increasing magnetic 
field strength as oc B~ A and area as oc A~ 2 . These ar¬ 
rays seem potentially advantageous in probing the axion 
“sweet spot” of 10 10 < f a < 10 12 GeV. These are very 
preliminary considerations, but appear worthy of further 
more detailed analysis. 

We note that our present paper and p] has not been 
without some controversy. In part, in some earlier ver¬ 
sions we mis-stated how the decoupling of the axion 
works for the OEDM in the Introduction. Hoswever, in 
[26^ a calculation was done in a pure Coulomb poten¬ 
tial, which yielded zero and was used to argue that our 
overall result is zero. As we have shown explicitly, the 
Coulomb potential part of the result is always a total 
divergence, and indeed vanishes. The OEDM is intrinsi¬ 
cally time dependent and is nonzero as an action. The 
axion decoupling has been thoroughly explored in the 
present paper and this issue has now been resolved fa¬ 
vorably for the present work and Q. (see HI). We have 
explicitly addressed the criticism that our result did not 
display decoupling, in the m a —> 0 limit. The decoupling 
is present, even in the radiation formula of eq. ©■ but it 
is not always manifest, and is subtle in a manner similar 
to the anomaly itself, as discussed above in Section III. 
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Appendix A: Slab Radiator Solution 

Here we give the solution for the electric dipole radia¬ 
tion form a conducting slab with an imbedded magnetic 
field. The power from a slab array can be estimated from 
the continuum limit of a discrete array of magnets. The 
emitted power from a discrete array of N magnets is of 
order: 


P = gl {HBohrf m A a elN (A.l) 

The axion Compton wavelength is m a ~ (27r/A). We can 
identify the magnetic field in the discrete array with 

(A.2) 


and the discrete array power is then: 

P ~ g\ (t^rn-XN 

~ g \ (A) #0 ( A / 27r ) 2 N ~ 9 a (A) 00 x ( area ) 

(A.3) 


Let us derive this result carefully from Maxwell’s equa¬ 
tions. The sourced Maxwell equation is: 

<9 t if r = gl^odtO + it (A.4) 

where ifo is a constant background magnetic field; 9 is 
oscillating axion field and if r and if r are oscillating re¬ 
sponse fields. 

We assume we have vacuum for x < 0 and a conducting 
medium, such as copper, for x > 0. B 0 and 9 perme¬ 
ate the vacuum and the conductor. The magnetic field, 
if 0 = Bqz is parallel to the plane of the conductor. We 
assume Ohm’s Law in the conductor 

~f = — crif (A.5) 

Introduce a vector potential in radiation gauge for the 


response 

fields: 




if r = 

-d t X 

B r — X x X 

X-X = 0 

(A.6) 

hence: 





% 

X -V : 

2 X = gX 0 d t 6 + X 

x > 0 

(A.7) 


d?X 

- v 2 ! = g X 0 d t e 

x < 0 

(A.8) 


The axion induced particular solution is: 

if 0 = —ffifo0 = — d t ~fo hence, <9 t 2 xto = glt 0 d t 9 

(A.9) 

Note x if r o = g x ifo = 0 since if 0 is assumed 
constant (or source free). Hence, V 2 yf 0 = 0 and ~fo is 
non-propagating and constant in space (curl-free). 

It is convenient to complexify the fields and identify 
the physical fields with the real parts. We write 9 => 
9q exp ( im a t ), then gltgdtO = im a g~$o9o exp ( im a t ), and 
ito = — igl^o9o/ma exp (im a t). The equation of motion 
becomes: 

df~A — V 2 ll = im a g~f o0oexp( im a t ) + ad fX {x > 0) 
d 2 X — \/ 2 X = im a g~X 0 9 0 exp (im a t) {x < 0) (A.10) 

X is composed of homogenous propagating and non¬ 
propagating solutions. Define: 

X = Xe irriat + Xe im “ t+k ' x - l9 ^° 6a e im at r x > 0 ) 

m a 

= J^ e im ^+ kx - 


0 rs “' f^Bohr (27r/A) ~ f^Bohr (j^a) 


m, 


{x < 0) (A.11) 
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We introduce a parameter: 7 = m a /<J and we will work We therefore have the resulting solution for small 7 << 1: 
in the limit of large conductance, 7 << 1. Substituting 
into the equations of motion we have: 


- (m 1 2 a + im 2 J 7 ) = gt Q 6 Q {m a /~i) 

-(ml + iml/'y-k ,2 )l$ = 0 

k 2 )~t = 0 


Hence: 


- \m n - . 


t = 


igl^odo 


m a (1 - * 7 ) 
k! = ±m a \/l + i /7 
k = m n . 


(A.12) 


(A.13) 


! 


1 


g B 0 6^jrn a t+i^x-^ 
9 - 3 q9q imat 


^ e i7r / 4 5 6 7 _|_ 


m a 


(x > 0) 


*3 




(l + v^e i7r/4 ) e (im “ (t+x)) 

*3^ofr 


£ e (im„t) (3 < 0 ) ( A . 16 ) 


The sign of £/ is determined by the asymptotic bound¬ 
ary condition that 7^ —> 0 as x — > 00 . The boundary 
conditions at x = 0 match the fields and their spatial 
derivatives on L and R: 

ig^oQo + ^ = ^ 

TO a (1 - *7) 

y/l + */ 7 ^ = ^ (A. 14) 

whence: 


We thus see, in the 7 —>■ 0 limit, the overall amplitude of 
—> 0 for x > 0 , and the propagating component atten¬ 
uates into the conductor with a skin depth of 1 / y/m a <x- 
For x < 0, in the 7 — >■ 0 limit, we see that the induced 
propagating plane wave has an overall amplitude that has 
become locked to the magnitude of the non-propagating 
axion induced particular solution. For x = 0 w see that 
the solution approaches the usual conducting boundary 
condition, ~3(x = 0 ) = 0 
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(A.15) 


The main point is to ilustrate that a conductor con¬ 
verts the non-propagating particular solution into de¬ 
tectable radiation. This is how the walls of the RF cavity 
pump radiation into the volume; the magnetic field in the 
empty volume plays no role, only its permeating the walls 
induces the signal. This also illustrates how the slab unit 
cell acts as a radiation source. 
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